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Abstract 

The stability of magnetized strange quark matter (MSQM) is studied in the MIT bag model 
with the density dependent bag pressure. In the consistent thermodynamic description of MSQM, 
the quark chemical potentials, the total thermodynamic potential and the anisotropic pressure 
acquire the corresponding additional term proportional to the density derivative of the bag pressure. 
The model parameter space is determined, for which MSQM is absolutely stable, i.e., its energy 
per baryon is less than that of the most stable 56 Fe nucleus under the zero external pressure 
and vanishing temperature. It is shown that there exists the magnetic field strength H umax at 
which the upper bound B ^ on the asymptotic bag pressure = B{qb 3> go) (Qo being the 
nuclear saturation density) from the absolute stability window vanishes. The value of this field, 
Humax ~(l-3) • 10 18 G, represents the upper bound on the magnetic field strength, which can be 
reached in a strongly magnetized strange quark star. It is clarified how the absolute stability 
window and upper bound on the magnetic field strength are affected by varying the parameters in 
the Gaussian parametrization for the density dependence of the bag pressure. 
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I. INTRODUCTION AND BASIC EQUATIONS 


After the conjecture that strange quark matter (SQM), composed of deconhned u, d and s 


quarks, can be the ground state of matter 


3j, it became the subject of intense researches. 


In the astrophysical context, this would mean that the formation of strange quark stars, 


made up entirely of SQM and self-bound by strong interactions, is possible 


4-6]. The birth 


of a strange quark star can proceed via conversion of a neutron star as a strong deflagration 
process during a few milliseconds Q, accompanied by a powerful neutrino signal 8]. If SQM 
is metastable at zero external pressure, it can be encountered in the cores of heavy neutron 
stars where the density of about several times nuclear saturation density can be sufficient 

n 

for the deconfinement phase transition to occur [9j. Such stars, composed of the quark core 
and hadronic crust, are called hybrid stars. Modern astrophysical observations, including 
data on the masses and radii, spin-down rates, cooling history, glitches and superbursts, do 
not disprove the existence of quark matter in compact stars. 

The other important feature is that compact stars can be endowed with the strong mag¬ 


netic held [l(J. Near the surface of magnetars - strong 


11 


y magnetized neutron stars - the 


held strength can reach values of about 10 14 -10 15 G 
up to 10 19 -10 2 ° G can potentially occur in the cores of neutron stars 


12]. Even stronger magnetic helds 


13]. The large pulsar 


kick velocities due to the asymmetric neutrino emission in direct Urea processes in the dense 


core o: 


helds 


a_ magnetized neutron star could be the possible imprint of such ultrastrong magnetic 


14 


17]. The origin of magnetar’s strong magnetic helds is yet under discussion, and, 


among other possibilities, it is not excluded that this can be due to spontaneous ordering of 


nucleon 


18 


19], or quark 20| spins in the dense interior of a neutron star. 


Strong magnetic helds can have significant impact on thermodynamic properties of cold 


dense matter 


3 - 29 ] 


In particular, the pressure anisotropy, exhibited in the difference 


between the pressures along anc 


strongly magnetized matter 


30 


perpendicular to the magnetic held, becomes relevant for 


35j. In this study, I consider strongly magnetized SQM 


(MSQM) taking into account the effects of the pressure anisotropy. I aim at finding the 
model parameter space for which MSQM is absolutely stable, i.e., its energy per baryon is 
less than that of the most stable nucleus 56 Fe under the zero external pressure and vanishing 
temperature. For the parameters from this absolute stability window, the formation of a 
strongly magnetized strange quark star is possible. 
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Note that in order to describe the confinement property of quantum chromodynamics, in 


the conventional MIT bag model 


this is achieved by introducing the density independent 


bag pressure by which quarks are confined in a finite region of space called a ’’bag”. The 
standard thermodynamic equations can be used to study quarks confined to a bag. Another 
phenomenological way to describe the quark confinement is to consider the density depen¬ 


dent quark masses 37-140]. In this case, an important issue of thermodynamic consistency 


arises. Because of the density dependence of the quark masses, the quark chemical poten¬ 
tials acquire an additional density dependent term and become effective 1411. In fact, such 


thermodynamically consistent approach was developed in the late work 


42], Note also that 


the quark confinement was modeled recently by the density- and isospin-dependent quark 


masses 


43 


27], 


hese phenomenological QCD models were applied to study MSQM in Refs 


29 


35 


in Refs. 121 



• 0 , 122 , 


24, 


45 J46L In particular, the effects of the pressure anisotropy were disregarded 


22], 24, [27j, in Refs. 45i, 146] only the matter contribution to the pressure 


anisotropy was considered, and both the matter and field contributions were accounted for 
in Refs. [35, 44], I 11 this study, I consider the absolute stability of MSQM in the MIT 
bag model with the density dependent bag pressure B(q b ). The advisability to extend the 
conventional MIT bag model came from the necessity to reconcile the different constraints 
on the bag pressure at low and high baryon densities obtained from heavy-ion experiments 
at CERN-SPS and astrophysical observations of neutron stars with the masses well above 
the mass of a canonical neutron star M ~ 1.4M 0 (M 0 being the solar mass) [47]. In other 
model frameworks, the density dependent bag pressure was used in Refs. |48|, 490 - Note that 
in the extended MIT bag model, the quark chemical potentials acquire the term proportional 
to the density derivative J^-, and the issue of the consistent thermodynamic description of 
MSQM becomes relevant. To that aim, I explore the formalism of the work 42], developed 
initially to describe the quark confinement by the density- (and/or temperature-) dependent 
quark masses, and, after proper modification, apply it to the case of the extended MIT bag 
model. 


Here MSQM will be considered as an uniform matter permeated by an external uniform 
magnetic field. I 11 the bag model, the matter part of the total energy density (excluding the 
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magnetic field energy contribution) reads 


E m = 0 


E 


Q'iQii 


(i) 


where 




( 2 ) 


and 0° is the thermodynamic potential for free relativistic fermions of ith species (i — 
u, d ;_s, e) in the external magnetic held, which is given by the same expression as in Refs. 


21 


122 , 


351 ] with the only difference that, according to the approach of Ref. 


42], the real (non- 


renormalized) chemical potentials /i,- should be substituted there by the effective (renormal¬ 
ized) chemical potentials /i,. 

At the given H , the differential form of Eq. (JTJ) is 


{dE m ) H (^^m) E ^ ^ h* d>Qi A ^ ^ Qi dfli, 


( 3 ) 


where 


With account of equation 


Eq. ([3]) acquires the form 


( dQ m) H = E ^ + E ^7 d ^ 

i ' i 

/da 0 ' 

Qi 


dfli j h' 


dB 


(dE m ) H — y A/h + J dg 


On the other hand, the fundamental thermodynamic relation at zero temperature reads 


( 4 ) 


( 5 ) 


50] 


( dE m ) ^ ^ l^idQi. 


( 6 ) 


By comparing Eqs. (J5]) and (J6]), and taking into account expression for the baryon number 
density qb = |(fti + £ , d + £?s), one gets the relationship between the real and effective chemical 
potentials 


_ - ^ 1 dB _ 

he he, h/ h/ E g ^ , J U, d, S. 


( 7 ) 
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Further I study charge neutral states of MSQM and assume that the chemical equilibrium 
with respect to weak processes is established among the fermion species with the corespond¬ 
ing conditions on the real chemical potentials [5, |42j. Note that, in view of Eqs. (I7|) . the 
effective chemical potentials fi, satisfy the same equations, as the real ones 


I hi Pu T pe~ j pd, /ts ■ 


( 8 ) 


The Hugenholtz - van Hove theorem establishes the thermodynamic relation between the 
pressure and energy density at zero temperature for nonmagnetized fermion matter 51]. For 


magnetized 
strength [30 


ermion matter, the total pressure is the anisotropic function of the magnetic held 


35]. In particular, the longitudinal p l and transverse p l pressures are differen 


By comparing expressions for the longitudinal pressure p l and energy density hq 
one can get the Hugenholtz - van Hove theorem for magnetized matter in the form 


35|, 


Pm — T ^ ^ Pi Qii (9) 

i 

where p l m is the matter part of the longitudinal pressure. I will preserve this equation also 
for MSQM in the extended MIT bag model. With account of Eq. ©. Eq. © takes the form 

_ QJB) 

Pm = ~^m T ^ 'X pi ~ Pi)Qi = ~^m T QB ■ (10) 

i 

At zero temperature, the matter part of the thermodynamic potential determined 
according to the standard thermodynamic equation 


E m y ] Pi Qj , 

i 

with account of Eq. (JTj) becomes 

= — y ^ (pi — Pi) Qi = Qn 


( 11 ) 


Qb 


dB 

dpB 


( 12 ) 


By comparing Eqs. m and (Eli. one arrives at the thermodynamic relationship p l m = 


— Vt m . The mat 
equation 


31, 


44 


;er parts of the longitudinal and transverse pressures are related by the 
52] 

Pm Pm = HM : (13) 


where M — —eetf 1 is the system magnetization. After summarizing Eqs. (HJ, I©, l l If II I, 
03 , and accounting for the pure magnetic held contribution, the total energy density E, 
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longitudinal p l and transverse p l pressures for MSQM with the density dependent bag 
pressure can be written in the form 


H 2 


E — + piQi) + — + B , 


l \ ^ nO 95 „ <95 


pl = _ En o-™ + |!. B+ ^ 


(14) 

(15) 

(16) 


In the case of the density independent bag pressure, Eqs. (1X41) — (jT6|) go over to the cor¬ 
responding equations of Ref. J 35 I - Because of the breaking of the rotational symmetry by 
the magnetic field, the longitudinal p l and transverse p t pressures are not the same. There 
are two different contributions to the pressure anisotropy: the matter contribution propor¬ 
tional to the magnetization M, and the magnetic field contribution given by the Maxwell 


term El. Note that it was argued recently 


87 r 


53| (and discussed before in Refs. [54|, [55]) that 


the magnetization contribution to the energy-momentum tensor is canceled by the Lorentz 
force associated with the magnetization current, and, taking that into account, the pres¬ 
sure anisotropy doesn’t occur. In the given case, the direct answer on this question is that 
I consider a spatially uniform distribution of the magnetic field and matter density, and, 
hence, the magnetization current density jj = V x M (in units with c = 1) is exactly zero. 
Therefore, the associated Lorentz force density j& x H = 0, and cannot compensate the term 
with the magnetization in the transverse pressure. The other point is that the main, most 
principal source of the pressure anisotropy is provided by the pure magnetic field contribu¬ 
tion 5 2 /87 t, and there is no compensating effect for it. Just this term, as will be shown 
later, plays the main role in establishing the upper bound for the magnetic field strength in 
a magnetized medium. 

For a spatially nonuniform case, the argument, raised in 531] . is based on the consider¬ 
ation of the balance of the volume forces in the stationary state: <9 {T kl = 0, T kl being the 
spatial components of the quantum-statistical average of the energy-momentum tensor of 
the system. Under writing the condition <9 \T kl = 0 with account of the Maxwell equation for 
the electromagnetic field tensor, the Lorentz force density associated with the magnetization 
current density, cancels the term with the magnetization, and the corresponding first inte¬ 
gral doesn’t contain the magnetization as well. While this step is correct, it is then wrongly 


concluded in Ref. 


53] that the system’s equilibrium (stationary state) is determined by the 
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thermodynamic pressure and any anisotropy in the matter pressure does not appear. Note 
that the pressures in the system are determined as the spatial diagonal elements T kk of the 
energy-momentum tensor in the system’s rest frame. In the stable stationary state, they 
should be positive T kk > 0. The conditions T kk > 0 and diT kl = 0 are, obviously, different 
and both should be satisfied in the stable stationary state. 

In fact, the consistent derivation shows that the above mentioned cancelation occurs, 
with account of the Maxwell equation, only in the derivative diT kl , but not in the energy- 
momentum tensor itself. An example of the spatially uniform case, discussed above, clearly 
confirms that. The introduction by hand of the Lorentz force contribution associated with 
the magnetization current to the energy-momentum tensor is the artificial step, not con¬ 
firmed by the consistent derivation. Just the opposite, the consistent microscopic deriva¬ 
tions |3ll . |52j show that such a contribution is missing in the energy-momentum tensor, 
and, hence, the anisotropy in the matter part of the total pressure is present as well as the 
anisotropy caused by the contribution of the magnetic held. 


II. NUMERICAL RESULTS AND DISCUSSION 

Now I will determine the absolute stability window of MSQM, subject to charge neutrality 
and chemical equilibrium conditions, at zero temperature. The equilibrium conditions for 
MSQM require vanishing the longitudinal p l and transverse p l pressures 

dQ 

^ = -0 = 0 , p t = -n + H— = 0, (17) 

where fl = VL m + is the total thermodynamic potential of the system. In order to be 
absolutely stable, the energy per baryon of MSQM should be less than that of the most 
stable 56 Fe nucleus under the equilibrium conditions (fTTj) . On the other hand, at H = 0, 
the experimental observation proves that two-flavor quark matter, consisting of u and d 
quarks, is less stable compared to 56 Fe nucleus at zero external pressure and vanishing 
temperature js|. I will also retain this constraint for strong magnetic fields H > 10 17 G, 
although, strictly speaking, it is unknown from the experimental point of view whether two- 
flavor quark matter is less stable than 56 Fe nucleus under the equilibrium conditions (TT7T) in 
such strong fields. Thus, for determining the absolute stability window of MSQM, I will use 
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the following constraints: 


J - J m 

Qb 


< e^( 56 Fe) < ^ 

uds 


(18) 


ud 


Regarding, for the rough estimate, 56 Fe nucleus as a system of noninteracting nucleons, 
magnetic fields H > 10 20 G are necessary in order to significantly alter its energy per 
nucleon e#( 56 Fe) 23]. Since 1 will consider magnetic fields H < 5 ■ 10 18 G, further I use the 
approximation e_H-( 56 Fe) m e 0 ( 56 Fe) = 930 MeV. 

For determining the absolute stability window from constraints (TT7T) and (fT8|) . 1 utilize 


the Gaussian parametrization for the bag pressure as a function of the baryon density 

B(q b ) = Boo + (B 0 - Roo) e ~ P , 


m- 


(19) 


where B 0 = B(qb = 0), B^ = B(qb Qo) (go — 0.17 frn -3 ), and the parameter (3 controls 
the rate of decrease of the bag pressure from B 0 to B^. In numerical analysis, I use [3 = 0.01 
and f3 = 0.17. Note that the equilibrium conditions (1T7T) contain the derivatives |j|, g ^g H ■ 
Further I will assume that the magnetic field affects the bag pressures B o and B^ in the 
same way, and, hence, the difference AB 0 = B 0 — B^ is independent of H and will be 
considered as a parameter. Then the actual dependence of the bag pressure B on H is 
contained in the asymptotic bag pressure B <*,: |j| = while d ^ gH = 0. The absolute 
stability window will be built in the plane (H, B 0Q ) under various fixed values of ABq and 
f3. The upper (lower) bound B ^ (B 1 ^) on the asymptotic bag pressure B 0Q can be found 
from the first of the equilibrium conditions (Tl7T) : 


b u j \ h ) = - Y. ^-fr (20) 

i=u,d,s,e 

(i=u,d,e) 



after finding the effective chemical potentials /b from the first constraint to the left (right) 
in (TT8|) . taken with the equality sign, and charge neutrality and chemical equilibrium con¬ 
ditions (J8]). In Eq. (l20jh the baryon density Qb should be determined from Eq. (J4J) after 
finding the chemical potentials /b at the given H. 

Fig. 1 shows the dependences B^(H) and B^H) for the current quark masses rn u = 
rrid = 5 MeV and m s = 150 MeV. The upper bound B ^ stays, at first, practically constant 
and then, beginning from the magnetic field strength H somewhat smaller than 10 18 G, 











FIG. 1. (Color online) The absolute stability window in the plane (H,B 00 ) for MSQM at zero 
temperature with B(qb) given by Eq. (fl9l) at (a) f3 = 0.01 and (b) /3 = 0.17, and with a variable 
parameter ABq. The upper B ^ and lower B ^ bounds are shown as the upper and lower curves, 
respectively, in the pairs of the similar curves (see also comments in the text). 

decreases. For example, the maximum value of B corresponding to H = 0 (which is prac¬ 
tically indistinguishable from the value of B ^ at H = 10 16 G) is B ^ max ps 74.9 MeV/fm 3 for 
AB 0 = 0, independently of the value of /3; for A B 0 = 75 MeV/fm 3 , B^ max ps 3.2 MeV/fm 3 
at (3 = 0.01, and ~ 49.1 MeV/fm 3 at (3 = 0.17. The upper bound B ^ vanishes 

at H umax ps 3.1 x 10 18 G for A B 0 = 0 at any /3; for A B 0 = 75 MeV/fm 3 , B ^ vanishes 

at H umax ps 1.1 x 10 18 G for [3 = 0.01, and at H umax ps 3 x 10 18 G for (3 = 0.17. In 

stronger magnetic fields H > H umax , in order to satisfy the equilibrium conditions, the 
upper bound B ^ on the asymptotic bag pressure B x had to become negative, contrary to 
the constraint B^ > 0. This means that, under the equilibrium conditions and in magnetic 
fields H > H umax , MSQM cannot be absolutely stable. 

The behavior of the lower bound B^^H) is similar to that of B^^H). At H = 0, the 

maximum value of B 1 ^ (which almost coincides with the value of B^ at H = 10 16 G) 

is B l 0o max ps 56.5 MeV/fm 3 for A B 0 = 0, independently of the value of f3; for A B 0 = 
25 MeV/fm 3 , B' oo max ps 32.0 MeV/fm 3 at (3 = 0.01, and max ps 40.0 MeV/fm 3 at 
f3 = 0.17. The lower bound B^ stays practically constant till magnetic helds somewhat 
smaller than 10 18 G, beyond which B^ decreases. The lower bound B^ vanishes at Hio ps 
2.7xl0 18 G for A B 0 = 0 at any f3; for AB 0 = 25 MeV/fm 3 , B ^ vanishes at H i0 ps 2.2x10 18 G 
for f3 = 0.01, and at H[ 0 ps 2.5 x 10 18 G for j3 = 0.17. Under the equilibrium conditions and in 
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FIG. 2. (Color online) The dependences Qb(H) for MSQM and magnetized two-flavor quark matter 
(the upper and lower curves, respectively, in the pairs of the similar curves), determined under the 
equilibrium conditions at E m /gi = 930 MeV for the same values of the parameters ABq and /3 as 
in Fig. 1. The bounding vertical lines on the right correspond to H = H umax . 

the fields H > Hi 0 , the lower bound B^ would be negative. Because B^ > 0, the inequality 
5oo > B ^ would be fulfilled always in the fields H > Hiq. Thus, in order for MSQM to be 
absolutely stable, the magnetic held strength should satisfy the constraint H < H umax . In 
fact, the value H umax represents the upper bound on the magnetic held strength which can 
be reached in a magnetized strange quark star. Note that the upper bound B ^ decreases 
with increasing the parameter ABq at the given j3. Hence, B ^ vanishes at a smaller magnetic 
held for a larger value of A B 0 , i.e., the upper bound on the magnetic held strength H umax 
decreases with increasing A B 0 at the given j3. On the other hand, the upper bound B ^ 
increases with increasing the parameter (5 at the given ABq. Hence, the larger /3 is, the 
larger the upper bound H umax at the given A B 0 . 

Note also that the lower bound B^ decreases with increasing the parameter A B 0 . Above 
certain value of A B 0 , which depends on f3, the lower bound B^ would become negative, 
and, hence, for B^ > 0 the inequality B > B^ would be always fulfilled at any H. For 
such ABq and (3, the absolute stability window corresponds to 0 < B^ < B This is just 
the case for A B 0 = 75 MeV/frn 3 at j3 = 0.01, and for A B 0 = 125 MeV/fm 3 at j3 — 0.17, 
shown in Fig. la and Fig. lb, respectively. 

Fig. 2 shows the dependences Qb(H) for MSQM and magnetized two-flavor quark matter, 
determined under the respective equilibrium conditions at E m j = 930 MeV, for magnetic 
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fields 77 < H umax . In fact, the corresponding lines g B (H) and g l B (H ) represent the upper 
and lower bounds on the baryon density to ensure the absolute stability of MSQM. It is 
seen that the upper and lower g l B bounds stay practically constant till the magnetic 
held strength being somewhat smaller than 10 18 G, and then increase till it reaches the 
corresponding maximum value. The increase of the parameter ABq at the given /3, and the 
increase of the parameter f3 at the given A7> 0 lead to the increase of both the upper g B and 
lower g l B bounds. For example, at H = 0 the minimum value of g B (being almost the same 
as the value of g B at H — 10 16 G) is g Bmin ~ 2.1 g 0 for AB 0 = 0, independently of f3; for 
AB 0 = 75 MeV/fm 3 , g Bmin ~ 2.2 g 0 at (3 = 0.01, and g Bmin ~ 2.9po at (3 = 0.17. Similarly, 
at H = 0 the minimum value of g l B is g l Bmin ~ 1.4 g 0 for AB 0 = 0, independently of /3; for 
AB 0 = 75 MeV/fm 3 , g l Bmin ~ 1.5p 0 at (3 = 0.01, and g l Bmin ~ 2.2p 0 at /3 = 0.17. Note 
that magnetic helds 77 ~> 10 18 G strongly affect the upper and lower g l B bounds from the 
absolute stability window. 

In conclusion, I have considered MSQM under charge neutrality and chemical equilib¬ 
rium conditions in the MIT bag model with the density dependent bag pressure B(g B ). I 
aimed to determine the range for the magnetic held strength 77, asymptotic bag pressure 
Boo = B(g B p 0 ), and baryon density g B , for which MSQM is absolutely stable, i.e., its 
energy per baryon is less than that of the most stable 56 Fe nucleus under the zero pressure 
conditions (fT71) and vanishing temperature. In fact, this requirement sets the upper bound 
on the parameters from the absolute stability window. The lower bound is determined from 
the constraint that magnetized two-flavor quark matter under equilibrium conditions (TT7j) 
and zero temperature should be less stable than the most stable 56 Fe nucleus. This con¬ 
straint is extended from the weak terrestrial magnetic helds, where it has direct experimental 
confirmation, to possible strong magnetar interior magnetic helds 77 > 10 1 ' G, where such 
confirmation is wanting. It has been shown that there exists the magnetic held strength 
Humax at which the upper bound B ^ on the asymptotic bag pressure Boo from the absolute 
stability window vanishes. In fact, the value of this held, H umax ~(l-3) • 10 18 G, repre¬ 
sents the upper bound on the magnetic held strength, which can be reached in a strongly 
magnetized strange quark star. I have studied the effect of the parameters in the Gaussian 
parametrization for the bag pressure on the absolute stability window and upper bound 

Hy m ax • 

It is interesting to note that the obtained estimate for H umax in strange quark stars is 
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similar to the estimate H ~ (1-3) • 10 18 G for the maximum average magnetic fields in stable 


neutron stars, composed of strange baryonic matter 


bound H umax in strange quark stars may be further improved 
bag model the effects of the perturbative quark interactions 


56]. The found estimate of the upper 
Dy including within the MIT 


57]. It would be of interest 


also to extend this research to the case of the spatially nonuniform magnetic held, whose 


realistic proh 


equations 


58 


e should be determined from the solution of the coupled Einstein and Maxwell 
. The other interesting problem is to take into account the vacuum corrections 


due to the magnetic held, which are, however, beyond the scope of the MIT bag model. To 


that aim, one can utilize the quark chiral models 
Ref. 


59h 61|. Nevertheless, as was shown 


m 


62], these corrections become noticeable only in strong magnetic helds H > 3 x 10 19 G, 


and, therefore, one can expect that their effect on the maximum magnetic held in strange 
quark stars will be of less importance. 
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